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RELATIVE-HUMIDITY  PARAMETERIZATION  OF 
THE  NAVY  AEROSOL  MODEL  (NAM) 


INTRODUCTION 

The  validity  of  the  relative-humidity  parameterization  of  the  Navy  aerosol  model  (NAM)  (1]  is 
evaluated  and  improved  in  this  report  for  several  reasons: 

•  The  three  lognormal  size  distributions  corresponding  to  the  three  components  of  NAM  consist 
in  part  of  particles  with  sizes  smaller  than  0.1  This  may  be  a  source  of  error,  because  the 
existing  humidity  parameterization,  based  on  a  relationship  given  in  Ref  2,  only  gives  accurate 
results  for  particles  greater  than  0.1  ^m.  Of  the  three  components  of  NAM,  the  continental  or 
gas-to-particle  component  is  most  susceptible  to  this  error,  since  most  of  its  particles  are  smaller 
than  0.1  /*m. 

•  The  accuracy  of  the  humidity  relationship  is  limited  to  relative  humidities  (RKs)  less  than 
about  98%  (21.  It  would  be  desirable  to  extend  the  validity  of  NAM  beyond  98%.  While  the 
RH  range  between  98%  and  100%  is  narrow,  and  thus  infrequently  found  in  the  maritime  atmo¬ 
sphere,  it  is  nevertheless  important,  because  of  the  large  changes  in  particle  size  at  those  hu 
midities. 

•  A  further  improvement  is  to  provide  a  separate  humidity  relationship  for  each  component  for 
which  a  different  particle  type  is  considered;  presently  Fitzgerald’s  relationship  (2)  for  clean 
maritime  aerosols  (sea  salt)  is  used  for  all  three  components. 

Fitzgerald’s  relationship  (2]  gives  the  dependence  of  particle  size  on  the  dry  size  of  the  particle 
and  on  RH.  It  is  an  approximate  form  of  the  exact  equation  13-5]  for  the  equilibrium  size  of  solution 
droplets.  In  arriving  at  his  relationship,  Fitzgerald  assumes  that  the  curvature  effect  (Kelvin  effect)  of 
the  particles  is  negligible.  This  results  in  a  numerically  convenient  formula  but  also  causes  the  relation¬ 
ship  to  lose  accuracy  outside  of  the  particle  size  and  humidity  ranges  previously  mentioned.  Other 
approximate  solutions  to  the  exact  growth  equation  have  been  proposed.  Included  are  those  found  in 
Refs  4,  5,  6,  7,  and  8.  Those  solutions  are  also  unable  to  accurately  predict  the  dependence  of  particle 
size  on  humidity  for  small  particles  and  for  high  RH.  Reference  5  gives  predictions  near  RH  ■■  100%.; 
however,  none  gives  an  analytical  expression  which  is  valid  for  all  the  particle  sizes  found  in  the  NAM 
size  distributions,  and  for  the  entire  interval  from  0%  to  100%. 

The  use  of  the  exact  growth  equations  for  a  purpose  such  as  NAM  is  also  undesirable  because  the 
equations  are  implicit  in  particle  size  so  they  cannot  be  solved  directly  for  particle  size  as  a  function  of 
RH.  The  equations  can  be  inverted  to  obtain  particle  size  with  an  error  minimization  technique  [9],  or 
tables  of  particle-size  values  can  be  prepared.  Neither  method  is  useful  for  NAM  where  a  simple  ex¬ 
plicit  relationship  which  uses  a  minimum  of  computer  time  is  preferred. 

An  improved  approximation  of  the  exact  growth  equations  is  described  here.  By  including  the 
influence  of  the  Kelvin  effect,  this  approximation  avoids  the  shortcomings  of  the  previous  approxima¬ 
tions.  It  gives  acceptable  accuracy  in  predicting  particle  size  as  a  function  of  RH  for  dry  particles  rang¬ 
ing  from  0.01  to  10  urn,  and  for  an  RH  range  including  all  humidities  up  to  100%.  The  particles  are 
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assumed  10  be  chemically  homogeneous,  the  hysteresis  effect  is  not  included,  and  the  results  are  given 
for  increasing  humidity. 

The  new  approximation  is  applied  to  the  lognormal  size  distributions  in  NAM  to  discover  their 
behavior  as  a  function  of  RH,  especially  for  the  high  values  of  RH  and  for  the  NAM  component  with 
the  large  fraction  of  particles  less  than  0.1  fun.  This  behavior  is  used  to  parameterize  the  geometric 
mean  and  standard  deviation  of  the  lognormal  distributions  in  terms  of  RH.  The  RH  dependence  of 
the  aerosol  optics  (volume  extinction  coefficient)  is  estimated  for  the  lognormal  size  distributions  of 
NAM  in  which  the  full  humidity  parameterization  as  well  as  simplified  versions  which  neglect  the  Kel¬ 
vin  effect  are  used.  Limits  are  given  for  the  applicability  of  the  humidity  parameterization  for  sea-salt, 
urban,  and  rural  aerosols;  and  conclusions  are  drawn  as  to  the  validity  of  Fitzgerald’s  relationship  (2). 

LOGNORMAL  SIZE  DISTRIBUTIONS 

The  form  of  the  lognormal  distribution  function  used  by  Gathman  (1)  for  NAM  is 

-  V  ~  exp  -C  f  In  r  -  In  (/  r0l)  1  3  ,  (1) 

dr  ft  f  l  J 

where  dN/dr  is  the  number  of  particles  per  cm3  per  micrometer;  /  gives  the  number  of  the  component 

in  NAM  with  the  size  of  the  number  in  the  same  sequence  as  the  size  of  the  geometric  means;  r  is  par¬ 

ticle  radius  in  micrometers;  r0/  is  the  mode  radius;  C  ”  1/(2  x  s2)  where  s  is  the  standard  deviation  of 
the  size  distributions;  A,  are  parameters  fitted  to  the  meteorological  conditions  and  particle  concentra¬ 
tions  of  the  NAM  components  (1);  and  /  is  the  RH  effect  given  by  the  ratio  of  the  particle  radius  at 
the  ambient  RH  to  the  particle  radius  at  the  standard  RH  of  80%. 


/ 


o  l'/3 
6(1  -  S)  ' 


(2) 


where  S  is  the  saturation  ratio. 


Rather  than  use  Eq.  (1)  in  this  study,  the  normalized  form  of  the  lognormal  distribution  function 
[10]  is  used.  It  is  given  by 

dN  3  N,  |  _  [logr  -  log  if  rtt)Y  1 

dlogr  (?n)i/2  logb,  CXP  \  2 Uogb,)2  }’ 

where  N,  is  the  number  of  particles  per  cm3  for  each  component,  b,  is  the  geometric  standard  devia¬ 
tion,  and  rgl  are  the  geometric  means  of  the  size  distributions.  Equation  (3)  was  used,  because  for 
.V,  —  1  integration  over  all  particle  sizes  gives  dN/ dlogr  —  1  which  makes  the  comparison  of  distribu¬ 
tions  calculated  for  various  values  of  RH,  bg ,  and  rgl  easier,  and  permits  the  comparison  of  the  present 
results  with  similar  work  which  is  also  usually  described  in  terms  of  bg  and  rg%  e  g.,  [9]. 

Transformation  equations  [11]  can  be  used  to  equate  Eqs.  (1)  and  (3).  When  this  is  done,  bt/s 
«  2.87,  rgJra,  -  1.65,  NJA,  -  1.38  x  rg,  ;  and  bg  “  2.03,  rg]  -  0.05  ^ra,  rg2  -  0.4  ^m,  and  rg)  - 
3  3  p.m. 


AEROSOL  GROWTH  EQUATION 


The  ratio  of  the  water  vapor  pressure  p'r  ,  over  the  surface  of  an  aerosol  consisting  of  a  solution 
droplet  to  the  vapor  pressure  p„  over  a  plane  water  surface  is  given  by  the  usual  expression  [4] 


(4) 
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where 


<r'  is  the  surface  tension  of  solution  droplet; 
p  is  the  density  of  solution  droplet; 
pd  is  the  density  of  dry  hygroscopic  particle; 

Rv  -  4.615  x  106  erg  g~'K~'  is  the  specific  gas  constant  of  water  vapor; 
T  is  absolute  temperature, 

r  is  the  radius  of  solution  droplet, 

rd  is  the  radius  of  dry  hygroscopic  particle; 
i  is  the  van't  Hoff  factor; 

Mw  is  the  molecular  weight  of  water;  and 
M,  is  the  molecular  weight  of  dry  particle. 


Under  the  assumption  that  the  volumes  of  the  water  and  the  dry  particle  are  additive  in  the  solu¬ 
tion  droplet,  the  mass  mw  of  water  in  (he  solution  droplet  is  given  by 

mw  -  ~  (r3  p‘  -  r}  pd)  -  (r3  -  rd ),  (5) 

so  that 

P  d 

- ”  71 - u — • 

(r3  -  rj)pw 

Using  Hand's  [5]  expression  for  a-': 

a-'  -  y„ [ar,  (T0)  +  a(T0  -  T)  +  b  md!  mw 1,  (7) 


where  w(7V  ”  75.6  dynes/cm  is  the  surface  tension  of  pure  water  at  the  temperature  Ta  -  273.16  K, 
a  —  0.153  dynes  1/cm,  and  ya  the  effect  on  tr'  due  to  surface  effects  is  assumed  to  equal  1.00.  The 
assumptions  are  made  that  p'  -  Pw  in  the  exponential  term  in  Eq.  (4);  that  the  solution  droplet  is  in 
equilibrium  with  the  saturation  ratio  S  of  the  atmosphere  so  that  S  —  pr'l p\  that  the  linear  mass 
increase  coefficient  p  —  /  Mm/Ms  tabulated  for  the  sea-salt,  urban,  and  rural  particles  by  [5],  and  for 
the  ammonium  sulphate  particles  by  [12]  is  independent  of  the  rd;  and  that  the  dry  particles  are  chemi¬ 
cally  homogeneous. 


The  preceding  equations  are  combined  to  result  in  a  form  of  the  growth  equation 

exp((2cr')/(r  p».  Ry  D] 

1  +  m  mdf  mw 


(8) 


which  cannot  be  solved  directly  for  r  as  expected,  and  which  is  used  in  the  subsequent  numerical  calcu¬ 
lations. 


APPROXIMATE  GROWTH  EQUATION 

Values  of  5  were  calculated  by  introducing  a  large  number  of  closely  spaced  values  of  r,  one 
value  of  rd ,  constants  for  the  particular  aerosol  type,  and  values  ofp  into  Eq.  (8)  at  T  -  298  K.  Since 
Jl  is  itself  a  function  of  S,  an  iterative  procedure  was  used  to  choose  the  proper  values  of  p  in  applying 
Eq.  (8).  Figure  1  shows  the  results  of  some  of  those  calculations  for  ammonium  sulphate  solution  dro¬ 
plets.  Each  solid  line  gives  the  calculated  relationship  between  the  water  volume  mw  and  -log  S  for  a 
particular  ammonium  sulphate  rd.  By  neglecting  the  Kelvin  effect  the  relationship  between  mw  and 
-log  S  would  appear  in  the  form  of  straight  lines  in  Fig.  1,  which  is  approximately  the  case  for  large 
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Fig.  1  —  Volume  of  water  in  solution  droplets  containing  ammonium  sulphate  with  an  ini¬ 
tial  dry-particle  radius  as  given  in  equilibrium  with  the  atmospheric  water  vapor  saturation 
ratio  5.  The  exact  growth  equation  (solid  curves,  Eq.  (8)1  is  compered  with  the  approxi¬ 
mate  equation  (x.  Eq.  (1$)]. 


values  of  r  as  expected.  The  increasing  nonlinearity  of  the  curves  for  decreasing  values  of  rd  demon¬ 
strates  the  importance  of  including  the  Kelvin  effect  for  those  dry  particle  sizes. 

The  analytical  expression 


y 


c 


(9) 


(where  x  —  mw  4w/3 (r3  -  rd)  and  y  —  -log  S  )  was  Fitted  to  the  growth  curves  in  Fig.  1. 
Equation  (9)  is  explicit  in  r  so  that  it  can  be  solved  for  r  which  gives 

A  C  D 


r  — 


B  C  -  log  S 


+  ri 


1/3 


(10) 


where  D  *  3/4sr  . 

Figure  2  shows  that  power  laws  relate  A  and  C  to  rd  with  good  accuracy.  Straight  lines  fitted 
through  13  geometrically  equal  values  of  rd  between  0.01  and  10  fi m  give 


A  -  1.490  x  107r/510  , 

(11) 

C  -  1.352  x  10'7  rd  1  428  . 

(12) 

The  value  of  B  is  independent  of  rd  : 

B  -  2.30  x  10-4 

(13) 

D  -  0.2387. 

(14) 
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Fig  2  —  Relationships  or  A  and  C  to  the  dry  particle  radius  of 
ammonium  sulphate  particles  [Eqs,  (II).  (12)1 


Combining  Eqs.  (9M14)  gives  ihe  final  form  of  the  approximate  growth  equation  for  ammonium 
sulphate  particles: 


r 


Cl  r. 


C2 


C 3  -log  S 


(15) 


where  the  constants  Cl,  C 2,  C3,  C4  are  listed  in  Table  !.  A  similar  fitting  procedure  was  done  for 
the  sea-salt,  rural,  and  urban  aerosols;  Table  1  also  gives  the  appropriate  constants. 


Table  1  —  Constants  for  the  Growth  Equation,  Eq.  (15);  and  for  the 
Lognormal  Size  Distribution,  Eqs.  (19),  (22),  (27) 


Aerosol 

Model 

NAM 

Component 

rgd 

(cm) 

Cl 

C2 

C3 

C4 

C5 

C6 

Sea  Salt 

2 

2.460E-5 

2.0091 

0.7674 

3.079 

2  572E-I1 

-1.424 

10.23 

Sea  Salt 

3 

1.790E-4 

2.009) 

Urban 

1 

3.89 1 E — 6 

1,997 

0.3926 

3.101 

4.190E-II 

Hli[l 

IBS 

9.742 

Rural 

1 

4.210E-6 

2.002 

0.2789 

3.115 

5.4I5E— 11 

BE  :m 

urn 

9.681 

(NH4)2SO< 

0.4809 

3.082 

3  110E-11 

mEEM 

The  temperature  dependence  of  Eq.  (15)  is  found  primarily  in  C3  due  to  the  sensitivity  of  the 
Kelvin  effect  to  temperature;  C3  can  be  temperature  corrected  with  the  expression  C3(D  =*C3[1  + 
.004(298  -  T)  ). 

Figure  1  illustrates  the  behavior  of  Eq.  (15)  for  ammonium  sulphate  particles  by  values  of  mK  and 
-logS(x)  calculated  for  each  rd  curve.  There  is  agreement  between  the  exact  (Eq.  (8)1  and  the 
approximate  (Eq.  (15)]  growth  curves.  In  this  case  as  well  as  for  the  other  aerosol  types,  the  largest 
differences  occur  for  small  values  of  5  which  fall  below  the  deliquescent  point  of  the  particles.  In  this 
range  rapid  changes  occur  in  the  value  of  Jl.  The  mean  relative  error  in  %.  (  lOOx/lASlr (approx.)  -  r 
(exact)]/.-  (exact)],  between  the  exact  and  approximate  values  of  r  for  S  >  “0  75  and  for  each  rd 
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curve  was  calculated  for  all  aerosol  types  and  is  given  in  Pig.  3.  In  all  cases,  the  mean  error  is  less  than 
about  6%  for  the  range  of  rd  between  0.01  and  10  nm.  The  maximum  error  of  13.5%  was  found  for 
sea-salt  particles  with  rd  -  0.01  /im  at  S  —  0.75. 


Fig.  3  —  Mean  error  in  the  approximate  growth  equation 
ui  c„.ni  jred  10  the  exae:  growth  equation  as  a  function 
of  the  dry  particle  radius  and  the  different  aerosol 
models 


RELATIVE-HUMIDITY  PARAMETERIZATION  OF  THE  LOGNORMAL  SIZE 
DISTRIBUTIONS  OF  NAM 

Using  the  approximate  growth  formula  (Eq.  (15)],  the  behavior  of  the  lognormal  size  distributions 
((Eq.  (3)1  of  NAM  can  be  investigated  as  a  function  of  5.  Given  that  Eq.  (3)  is  in  terms  of  rd  ,  and 
the  dry-particle  geometric  mean  rgd  and  standard  deviation  bgd,  and  that  /  — 1,  the  particle  size  distribu¬ 
tion  as  a  function  of  S  and  r  is  given  by 

dN  _  dN 
dlogrj  dlogrdJ  dlogr, 

where 

dlogr dl  -  [log(rdJ+\l rdJ_\)\l2  -  constant  , 
dlogrj  -  l/og(ry+|/r,_,)]/2  , 
and  j  is  the  index  for  the  particle  radius. 

Figure  4  shows  the  results  of  applying  Eq  (16)  to  the  midsized  component  of  NAM  consisting  of 
sea  salt  with  rgd  -  0.246  ^m  and  bgd  -  2.009  (these  values  are  both  smaller  than  the  values  of  rs  and 
bg  given  previously  which  were  for  S  —0.8).  Several  changes  are  evident  in  ’he  size  distributions  in 
Fig.  4  as  S  increases:  the  maximum  value  of  dN! dlogr  shifts  to  larger  particle  sizes  and  at  the  same 
lime  decreases  somewhat,  and  the  width  of  the  distributions  increases.  The  increase  in  width  means 
that  dlogr  is  no  longer  a  constant  like  dlog  ( rd )  but  changes  with  S  This  clearly  shows  that  the  mean 
as  well  as  the  standard  deviation  of  the  distributions  are  a  function  of  S. 


(16) 

(17) 

(18) 
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Fig.  4  —  Dependence  of  the  lognormal  midsized  sea-salt  component  of 
NAM  on  the  given  saturation  ratios 


A  large  number  of  calculations  were  required  to  obtain  the  curves  in  Fig.  4,  because  100  values  of 
rd  were  needed  for  each  decade  in  particle  size  to  achieve  acceptable  accuracy.  That  makes  the  use  of 
Eq.  (16)  less  than  desirable  for  an  application  such  as  NAM.  For  this  reason,  an  attempt  was  made  to 
describe  the  distributions  in  Fig.  4  with  a  more  practical  expression  given  by  a  form  of  the  lognormal 
distribution 


dN  2  _ _  _  [ logr  -  logrt(S)]2 

dlogr  -f  (27r)1/2  logb,(S )  6  P  2[/0£&f(S)]2 


(19) 


where  both  rg  and  bg  are  functions  of  S.  This  equation  should  be  used  in  NAM  instead  of  Eqs.  (1)  or 
(3)  when  particle-size  information  is  desired  at  all  humidities  including  100% 


It  is  also  apparent  in  Fig  4  that  some  distortion  results  in  the  originally  lognormal  distribution  for 
dry  particles  as  the  value  of  S  increases.  This  comes  about  because  the  stronger  influence  of  the  Kel¬ 
vin  effect  for  the  small  dry  particles  prevents  them  from  growing  as  large  as  the  larger  dry  particles. 
Thus  before  Eq.  (19)  can  be  considered  for  inclusion  in  NAM,  it  is  necessary  to  obtain  an  estimate  of 
how  badly  the  original  lognormal  distribution  is  distorted.  This  can  be  done  by  calculating  rt(S )  and 
bt(S)  for  the  actual  distributions  given  by  Eq.  (16)  and  comparing  those  distributions  with  Eq.  (19). 
The  expressions  for  r,(S)  and  bg(S )  are 


|l  log',} 

rt  (S)  -  10 


X  dlogr, 

dlogr.  > 


(20) 


and 


bg(S)  -  10 


dS 
1  dlogr j 


dlogr jUogfj  -  /ojTj ( 5 )  1 2 


dlogr. 
dlogr,  J 


(21) 
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Figure  5  shows  the  comparison  of  the  actual  distributions  to  the  lognormal  distributions  which  are 
in  terms  of  rg{S )  and  bg(S).  The  distortion  is  relatively  small  for  S  <  .99,  it  reaches  a  maximum  near 
S  —  .9999,  and  is  negligible  for  5  —1.00.  Under  maximum  distortion,  values  of  r  predicted  by  Eq. 
(19)  are  accurate  near  maximum  values  of  dN/dlogr ,  these  values  lose  accuracy  as  dN/dlogr  decreases, 
and  they  are  in  error  by  about  30%  for  dN/dlogr  reduced  three  orders  of  magnitude  from  its  peak 
value.  This  value  of  S  for  the  maximum  distortion  only  applies  to  the  midsized  sea-salt  component.  It 
shifts  to  a  value  of  S  -0.999  for  the  component  with  the  smallest  particle  sizes  and  to  S  —  0.99999  for 
the  component  with  the  largest  particle  sizes  In  all  cases  the  distortion  is  negligible  at  S  -  1.00. 
While  this  distortion  is  undesirable,  it  occurs  at  values  of  S  which  will  seldom  be  used  due  to  their 
large  sizes. 

Equations  (20)  and  (21)  are  also  too  computationally  involved  to  be  practical  with  NAM;  hence, 
approximate  forms  were  found.  The  approximate  expression  for  rg(S)  is  simply  given  by  Eq.  (15)  writ¬ 
ten  in  terms  of  the  dry  geometric  mean  radius  rgd  : 


rg(S) 


Cl  # 

C3  rf*  -log  5 


(22) 


Equation  (22)  predicts  the  value  of  rg(S)  with  an  accuracy  of  better  than  3.2%  over  the  entire 
range  of  S  and  for  all  three  aerosol  components. 


An  approximate  expression  for  Eq.  (21)  is  not  as  straightforward,  because  of  the  complex 
behavior  of  bg(S)  as  a  function  of  S  as  shown  in  Fig.  6.  A  reasonable  fit  to  the  solid  lines  given  by  Eq. 
(21)  is  the  tangent  function  in  the  form 

v  -  -  M  +  B  tan  (Or)].  (23) 


where 


ys)-*„  i 

6f(l.0)  -  bgd(dry)  2 


(24) 


y  -  log  (1  -  S),  (25) 

and 

A  —  C5  logrgd  +  C6.  (26) 

B  —  0  65  and  C  —  0.88  are  empirical  constants,  and  C 5  and  C6,  listed  in  Table  1,  are  constants  which 
depend  on  the  aerosol  component. 


Given  that  in  all  cases  6*0.0)  -  bgd  =  0  .85,  Eqs.  (23)  to  (26)  are  combined  to  give  the  desired 
approximate  relationship  for  bg(S): 


bg(S)  —  bgd  +  0  300  arctan 


C 5  logrgd  +  C6+  log  (1  -  S) 
0.650 


+  0.415. 


(27) 


Figure  6  shows  the  predictions  of  Eq.  (27)  for  the  sea-salt  component  which  was  assumed  to 
apply  to  all  three  lognormal  size  distributions  of  NAM  The  accuracy  of  Eq,  (27)  is  better  than  0.05  for 
sea  salt  as  well  as  for  the  other  two  components  which  give  similar  results.  Figure  6  also  demonstrates 
that  the  change  of  the  geometric  standard  deviation  is  limited  to  large  values  of  RH.  For  the  lognormal 
distribution  with  the  smallest  particle  sizes  the  geometric  standard  deviation  must  be  corrected  with  Eq. 
(27)  for  values  of  S  greater  than  about  0.90,  while  for  the  other  two  distributions  the  correction  occurs 
at  higher  values  of  S. 
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OPTICAL  PROPERTIES  COMPUTED  FROM  THE  LOGNORMAL  SIZE 
DISTRIBUTIONS  OF  NAM 


Until  now,  this  analysis  has  dealt  only  with  the  RH  parameterization  of  the  size  distributions  of 
NAM,  with  Eq.  (19)  giving  the  final  form.  Since  the  ultimate  objective  of  NAM  is  to  predict  the  optics 
of  the  aerosol  particles,  this  section  investigates  the  accuracy  of  the  aerosol  optics  calculated  for  Eqs. 
(3)  and  (19).  Equation  (3)  is  the  transformed  version  of  Gathman’s  [1]  equation  (Eq.  (1)}  used  in 
NAM.  Equation  (19)  should  be  more  accurate  for  calculating  the  optical  effects  than  Eq.  (3);  however, 
it  has  not  been  made  clear  what  the  accuracy  of  the  optics  are  using  Eq.  (3).  This  accuracy  should  be 
belter  than  the  limits  imposed  on  the  applicable  range  of  RH  [Eq.  (2))  by  Fitzgerald  [2],  because  the 
optical  properties  are  primarily  determined  by  particles  greater  than  0.1  nm  for  which  the  Kelvin  effect 
has  less  of  an  influence. 


To  determine  the  accuracy  of  Eqs.  (3)  and  (19)  for  calculating  the  aerosol  optics  a  complete  Mie 
calculation  is  not  done.  Instead  the  cross  section  for  extinction  of  the  aerosol  particles  is  assumed  to  be 
proportional  to  the  particles’  area.  This  assumption  gives  an  estimate  of  the  upper  limit  of  the  error  in 
using  those  equations  for  the  aerosol  optics,  and  it  permits  much  shorter  and  simpler  calculations.  The 
area  of  the  particles  described  by  the  size  distribution  in  Eq.  (16)  is  calculated  and  is  assigned  to  be 
exactly  proportional  to  the  extinction  coefficient  of  the  particles.  Areas  (extinction  coefficients)  calcu¬ 
lated  for  the  approximate  expressions  of  the  size  distributions  are  then  compared  to  the  exact  extinction 
coefficient. 


The  comparison  of  the  Eq.  (19)  (including  Eqs.  (22)  and  (27))  to  Eq.  (16)  as  a  function  of  RH 
and  the  NAM  components  gives  the  following  results:  the  discrepancy  between  the  extinction  coeffi¬ 
cients  determined  from  those  equations  is  small  except  where  the  distortion  of  the  size  distribution 
given  by  Eq.  (19)  is  large.  This  occurs  at  very  high  values  of  RH  as  discussed  before.  In  this  region 
errors  on  the  order  of  15%  arc  found.  For  even  larger  values  of  RH  the  errors  decrease,  and  disappear 
altogether  at  RH  «■  100%. 

Equation  (19)  is  next  simplified  by  removing  the  Kelvin  effect  [found  in  the  term  C3  rfd*  of  Eq. 
(22)1,  and  by  assuming  that  6,(5)  is  simply  a  constant  equal  to  2.03.  The  particles’  area  resulting  from 
this  modified  form  of  Eq.  (19)  is  again  compared  to  the  area  of  Eq.  (16)  with  the  following  results: 
differences  do  not  exceed  about  10%  for  RH  <  99%  for  the  urban  and  rural  aerosol  (component  1),  for 
RH  <  99.9%  for  the  sea-salt  particles  (component  2),  and  for  RH  <  99  99%  for  the  sea-salt  particles 
(component  3). 


If  the  additional  approximations 


and 


rfd2  =  constant  x  rt}d 


(28) 


are  made  to  Eq.  (22),  then 


-log  S  =1  -  5 


r,  (5) 


1  + 


)l/J 

’ 


(29) 

(30) 


which  is  equivalent  to  Eq.  (7)  in  Fitzgerald’s  [2]  report,  and  where  y  is  defined  by  Fitzgerald  (2)  as  the 
air-mass  characteristic.  For  the  sake  of  comparison  with  his  values  of  y,  the  values  of  y  found  here  are 
0.825,  0.965,  0.282,  and  0.173  for  sca-salt  particles  (components  2  and  3)  and  for  urban  and  rural  parti¬ 
cles.  These  values  agree  reasonably  well  with  Fitzgerald’s  (2)  values.  It  is  of  interest  to  note  that  the 
value  of  y  differs  for  the  two  sea-salt  components  This  demonstrates  that  in  the  growth  curve  (Eq. 
(30))  the  air-mass  characteristic  y  depends  not  only  on  the  origin  of  the  particles,  but  also  on  their  size 
distribution 
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For  RH 


80%  (S  -  0.80)  Eq.  (30)  becomes 
r,(0.80)  -  rgd 


1  + 


0.2 


(31) 


Dividing  Eq.  (30)  by  Eq. 


(31)  eliminates  the  dependence  on  rgd  and  gives 


,  '.<S) 

J  rt  (0.8) 


Cl  -  S 
C8(l  -  S) 


t/J 


(32) 


which  is  the  form  used  by  Gathman  (1)  in  NAM.  Table  2  lists  the  values  of  Cl,  C8,  and  the  range  of 
validity  of  Eq.  (32)  for  the  various  NAM  components.  The  values  of  C7  —  2  and  C8  -  6  used  by 
Gathman  111  for  all  NAM  components  are  close  to  the  present  values  for  the  sea-salt  particles  in  NAM 
component  3.  The  assumptions  in  Eqs.  (28)  and  (29)  add  only  a  small  error  to  Eq.  (32).  This  means 
that  extinction  coefficients  calculated  with  Eqs.  (19)  and  (32),  which  give  an  expression  equivalent  to 
Eq.  (3),  have  an  accuracy  of  better  than  about  10%  for  the  range  of  validity  (upper  limit)  shown  in 
Table  2.  Thus  under  most  situations  found  in  the  maritime  atmosphere,  Eq.  (3)  [with  /  given  by  Eq. 
(32)]  can  be  used  instead  of  Eq.  (19).  The  constants  in  Table  2  must  be  used  with  Eq.  (3),  and  the 
range  of  validity  should  not  be  exceeded,  since  large  errors  are  found  at  greater  values  of  RH. 


Table  2  —  Constants  for  Eq.  (32),  and  the  Validity  of 
Eq.  (3)  for  Calculating  Extinction  Coefficients. 


Aerosol 

Model 

NAM 

Component 

Cl 

C8 

Range  of 
Validity 

Sea  Salt 

i  2 

1.83 

5.13 

RH  <  99.9% 

Sea  Salt 

3 

1.97 

5.83 

RH  <  99.99% 

Urban 

1 

1.28 

2.41 

Rural 

1 

1.17 

1.87 

RH  <  99% 

SUMMARY  AND  RECOMMENDATIONS 

A  new  approximation  formula  for  the  equilibrium  size  of  hygroscopic  particles  as  a  function  of 
RH  and  particle  type  was  developed.  This  formula  includes  the  influence  of  the  Kelvin  effect  which 
gives  the  formula  a  range  of  applicability  much  larger  than  previous  approximation  formulas.  It  applies 
to  particles  with  radii  extending  from  0.01  to  10  /*m,  and  for  a  RH  range  from  0%  to  100%.  The  new 
formula  was  applied  to  the  size  distributions  of  NAM  to  discover  the  effect  of  having  neglected  the 
Kelvin  effect  previously,  and  the  behavior  of  the  distributions  at  high  values  of  RH.  A  new  version  of 
the  NAM  size  distributions  [Eq.  (19)1  was  developed  where  the  geometric  mean  and  standard  deviation 
were  both  parameterized  as  a  function  of  RH.  It  is  estimated  that  the  error  of  calculating  aerosol 
extinction  coefficients  with  Eq.  (19)  does  not  exceed  about  20%  over  the  entire  RH  range  from  near 
0%  to  100%,  and  is  significantly  less  for  most  values  of  RH,  when  all  sources  of  error  are  considered. 
Sources  of  error  are  found  in  the  approximation  formula  for  particle  growth  primarily  due  to  the  non¬ 
linear  behavior  of  the  mass  increase  coefficient,  in  the  approximations  used  for  the  geometric  mean  and 
standard  deviation  of  the  size  distributions,  and  in  the  distortion  of  the  lognormal  distributions  at  very 
high  values  of  RH.  This  accuracy  is  sufficient  for  NAM,  since  much  greater  uncertainties  are  caused  by 
the  mass  increase  coefficients  given  by  Hand  [5]  which  reflect  average  conditions  and  thus  do  not 
include  the  large  fluctuations  in  the  coefficients  which  could  occur  as  a  function  of  time. 

It  was  demonstrated  fo"'.  the  humidity  parameterization  in  Eq.  (19)  reduced  to  the  one  presented 
by  Fitzgerald  (21,  [Eq.  (30)1  and  used  by  Gathman  [1],  (Eq,  (32)1.  The  reduction  was  accomplished  by 
applying  a  series  of  approximations  which  included  neglecting  the  Kelvin  effect  and  assuming  that  the 
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geometric  standard  deviation  was  independent  of  RH.  The  error  in  using  Gathman  and  Fitzgerald’s 
humidity  parameterized  size  distribution  (Eqs.  (1)  and  (2))  for  calculating  extinction  coefficients  for 
NAM  also  does  not  exceed  20%,  except  for  values  of  RH  which  are  larger  than  those  listed  in  Table  2. 

Specific  recommendations  for  changes  to  NAM  in  LOWTRAN  are: 

•  For  both  sea-salt  components  no  changes  are  required.  The  difference  in  the  values  of  the 
coefficients  (C7  and  C8)  in  tL growth  equation  [Eq.  (32)1  determined  here  and  those  used  by 
Gathman  [1]  in  Eq.  (2)  are  small  and  cause  negligible  error  in  aerosol  optics  calculations. 

•  For  the  continental  or  gas-to-particle  component  of  NAM,  C7  and  C8  in  Table  2  should  be 
used  with  Eq.  (2)  (see  Eq.  (32)].  Optical  coefficients  should  be  calculated  for  the  urban  and 
rural  aerosols,  and  criteria  established  for  their  use  within  NAM. 

•  The  routine  use  of  NAM  in  LOWTRAN  is  limited  to  RH  <  99%,  because  over  this  RH  range 
Gathman’s  (1983)  formulation  (Eq.  (1)1  can  be  retained  in  NAM,  and  the  humidity  parameter¬ 
ization  gives  valid  results  for  all  components. 

•  The  description  of  NAM  in  LOWTRAN  should  reference  this  report  to  indicate  that  it  gives 
approximation  formulas  which  extend  the  RH  dependence  of  the  sea-salt,  urban,  and  rural  par¬ 
ticles  to  RH  -  100%.  This  extension  is  appropriate  for  special  studies,  such  as  for  example, 
particle  behavior  under  prefog  or  precloud  conditions. 
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